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Exact Phase Solutions of Nonlinear Oscillators on 
Two-dimensional Lattice 

Tsunehiro Yokoi, Hiroyasu Yamada, and Kazuhiro Nozaki 

Abstract 



We present various exact solutions of a discrete complex Ginzburg- 
Landau (CGL) equation on a plane lattice, which describe target patterns 
and spiral patterns and derive their stability criteria. We also obtain 
similar solutions to a system of van der Pol's oscillators. 

1 Introduction 

Near the onset of the Hopf bifurcation in continuous media, a complex ampli- 
tude of distributed oscillators is described by a normal form equation called the 
complex Ginzburg-Landau (CGL) equation. Some exact solutions of the 1+1 
, dimensional CGL equation are known as hole, solitary-wave and shock solutions 

■ £Q. In the 2+1 dimensional space, numerical studies showed existence of inter- 
< ^ > esting solutions describing target patterns and spiral patterns |2 [3] . However, 

explicit and exact solutions representing target patterns or spiral patterns have 
not been constructed yet. 

■ Here we consider a system of oscillators which undergo the Hopf bifurca- 
tion on a two dimensional lattice instead of continuous media. Each oscillator 
is governed by the normal form equation of the Hopf bifurcation 4 and cou- 
ples diffusively and dispersively with its nearest-neighbor oscillators so that the 
system becomes a spatially discretized version of the 2+1 dimensional CGL (dis- 
crete CGL) equation^. Such a discrete CGL equation has also been studied 

■ on a network in connection with In this paper, we show that the discrete 
CGL equation has exact solutions describing various phase patterns consisting 
of target patterns and spiral patterns. Each spiral pattern solution needs a 
defect in its origin of spiral, where oscillators are absent. This defect may corre- 
sponds to the singularity of the origin in a spiral solution in continuous media. 
In the other hand, target pattern solutions do not have any defects. We derive 
the stability condition of the exact solutions to modulational perturbations and 
conduct some numerical simulations to confirm the stability criterion. Similar 
phase pattern solutions are constructed to a system of van der Pol's oscillators. 
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2 Nonlinear Oscillators on Two-dimensional Lat- 
tice 

Let us consider the following system of normal form oscillators to the Hopf 
bifurcation on a two-dimensional lattice. 

^™ = -(^ + ifMi)\W n \ 2 W n + (y r + ivi)W n + (A r + iXi) J2(Wj - Wn) (1) 

where W n are complex amplitudes of oscillators and fx r , fa, v r , Vi, X r , Xi are 
real constants; fi r and v r are assumed to be positive (the subcritical bifurca- 
tion) and X^~n denotes summation over nearest-neighbor oscillators at site n. 
Through this summation term, each oscillator has diffusive (A r ) and dispersive 
(Xi) interactions with its nearest-neighbors. Here, we call Eq.QJ a discrete CGL 
equation. 

By changing variables W n — > W n exp[— ifit] and t — > fat, and replacing 
coefficients as p\ = fa/ fa, p 2 — u r /fa, p 3 — X r /fa and p^ — X^j fa, Eq.Q is 
transformed into 



' /Un = -(l + ipi)\W n \ 2 W n +p 2 W n + (p 3 + i Pi )Y,(W J -W n ), (2) 



dt 
or 

^ = -(1 + iPi)\W n \ 2 W n + (pa - 4( P3 + ip 4 ))W„ + (ps + ip 4 ) J] W,-, (3) 
where j»2 > 0. We use Eq.® or Eq.© as a basic system in this paper. 

3 Exact Solutions 

In this section, we derive some exact solutions of the discrete CGL equation © , 
which represent a target pattern, a spiral pattern and various mixtures of those 
patterns. A crucial step to find exact solutions is to assume 

4 

^2 Wj = (for all sites n). (4) 

Then, Eq. J3J) takes the same equation for all n: 

= -(1 + i P i)\W n \ 2 W n + (p 2 - 4(p 3 + iPi))W n . (5) 

We choose an equilibrium solution of Eq. (J5J . 

W n = ^M+^(«)) 5 (6) 
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where A, uj are real constants satisfying 

A 2 = P2-4p 3 , (7) 
uj = -pi(p 2 - 4p 3 ) - 4p 4 , (8) 

and a functional form of the phase cj)(n) with respect to n is determined later. 
Because oi A 2 > 0, we need 

P2 - 4p 3 > 0. (9) 

Since P2 > 0, the condition 10 is automatically satisfied for p$ < (negative 
diffusion or diffusionless case). The solution has the same amplitude A at 
all cites but its phase could vary from site to site. The variation of the phase 
must be determined so that Eq.|@} is satisfied. 

4 

e^') = (for all sites n). (10) 

We consider a case that there are only two values for 4>(n), that is, (j)(n) = 
or 7r. At every site, we set phase values of two oscillators among four nearest- 
neighbor oscillators as 4>(n) = 0, while the other two phases are chosen to be it. 
Then, the condition (|10fl is fulfilled at all sites, that is, 



E 



= (for all sites n). (11) 



This choice of phase values is quite simple but produces a variety of phase 
patterns such as stripe patterns and target patterns as exact solutions. We can 
construct not only a single target pattern but also a combined pattern of some 
target patterns as shown in Fig.l. Although we present only a combined pattern 
of four target patterns, combined patterns of more or less target patterns are 
also possible. We can construct various combined patterns of target patterns 
indefinitely if the system size becomes larger. 

When the lattice has a defect, that is, an oscillator is absent at a site, we 
should modify the above discussion for four oscillators around a defect. To four 
oscillators around a defect, the coupling terms in Eq.Q are changed to 

3 3 
(p3+ipA)^2{Wj-W n ) = (p 3 +ip4)(J2 W 3- 3W n) 

= (pa + *Pa) I (W n + J2 Wj) -4W r \ , (12) 

where n indicates one of the four oscillators and the summation term in Eq. @ 
becomes 

3 

(p3 + iPi){W n + ^ i W j ). (13) 
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Thus, Eg. (|ll)fl is modified for each n-th oscillator around a defect as 



J4>(n) 



0. 



(14) 



Comparing Ea. H14|) to Ea. (|10l) . the term exp(i<fi(n)) in (|14|) appears replacing 
the term which comes from an oscillator at a defect site. Therefore, for the n-th 
oscillator around a defect, the defect behaves as a mirror, that is, the defect 
acts as if it is an oscillator with the same phase as the n-th phase. Owing to 
this flexible property of a defect, we can construct spiral patterns with 2 k arms 
around k x k combined defects. For example, we obtain a spiral pattern with 
two arms around a single defect {k = 1) and a spiral pattern with four arms 
around four combined defects (k = 2) as shown in Fig. 2 and Fig. 4, respectively. 
Two separated defects produce two spiral patterns with the same or opposite 
rotation directions as shown in Fig. 3. We can also construct various patterns in 
which spirals and targets are concurrent and show one of the simplest examples 
in Fig. 2. On a large enough lattice, we could construct virtually an indefinite 
number of such spiral and target patterns. 





Figure 1: A target pattern. Black points denote (f>(n) = 0, while gray points 
denote </>(n) = 7r. 
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Figure 4: A four-arm spiral pattern and 
a two-arm spiral pattern 
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4 Stability of Solutions 

Let us investigate a linear stability of the exact solution: 

W°(t) = Ae l ^ t+ ^ n)) {4>{n) = or tt), 

4 

^ e ^0) = o ( for all site n y 

Substituting a perturbed solution W n = W® + SW n into Eq.JSJ), we have 

^-SW n = iujSW n - (l + i Pl )A 2 5W n 
at 

4 

-(1 + i Pl )A 2 e 2l ^ t+ ^5W* + (p 3 + ipi) §W r ( 15 ) 

As <^(n) = or tt, the coefficient of the third term in the RHS of Eq. l|15|) becomes 
— (l+ipi)A 2 e~xjp[i2u>t], which does not depend on a site number n. Therefore, all 
coefficients of the linearized equation (|15J) are independent of n, which greatly 
simplifies stability analysis. Changing the variable as 8W n — > SW n exp[— iuit], 
Ea. is transformed into 

d 4 

-SW n = -(1 + i Pl )A 2 (5W n + SW*) + ( P3 + i P4 ) 5W 3 . (16) 

We consider plane- wave perturbations 

SW n = a{t)e tts = a {t)e l{k * x+k y y) , (17) 

where x — (x,y) indicates a lattice point and k = (k Xl k y ) is a wave-number 
vector. Then, the summation term of EcL ljlfif) reads 

4 

SWj = ae ts (e ik * + e lk y + e - ik °> + e -' fc » ) 

= ae k x (2 cos k x + 2 cos k y ) 
= 26W n K, 

where K = cosk x + cosk y . Hence, Eq. l(T^|) reduces to 

j t 8W n = -(l + i Pl )A 2 (6W n + 6W:) + 2( P3 + ip 4 )K5W n 

= (2( P3 + ip 4 )K-(l + ip 1 )A 2 )8W n -(l + ip 1 )A 2 SW; i , (18) 

or in terms of the real and imaginary parts of SW n 

d_ ( &[SW n ] \_ f P3 K-A 2 -p 4 K \ ( K[5W n ] \ 

dt I $*[SW n \ - Z \ p 4 K~ Pl A 2 p 3 K I I $$[SW n ] J ■ Uyj 
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Eigenvalues of Eq. (|19[) are 

2p 3 K -A 2 ± ^-Ap 2 K 2 + 4 Pl p 4 A 2 K + A 4 . (20) 
We denote an eigenvalue as a+ which has a larger real part than the other. 

a+ = 2p 3 K - A 2 + 9? pJ ' -Ap\K 2 + 4p lP4 A 2 K + A 4 . (21) 

If a+ < for any K in \K\ < 2, the various phase solutions are stable to a 
perturbation l|17|l . For stability, a + must take local maximum at K=0 since 
a + = at K = (see Fig|31 FigHJJ). The derivative of a + at K = is given as 



daj 



dK 



= 2 



K=0 



P3 ~Pi- 



2p 4 K -piA 2 



y/ -4p 2 K 2 + 4 PlP4 K + A 4 

2(j»3 +PlP4)- 



J K=0 



A condition p 3 + piPi = is necessary for stability. When p 3 + pip^ — 0, we 
have 



a+ 



K = 



= u 



3? 



= 0, 



^ +2 



PiPi 

P3 



A 2 



- -1 



A 2 



which shows that K = A 2 /2p 3 is also a vanishing point of a+, where its deriva- 
tive with respect to K is finite. For stability, we need 

\A 2 /2p 3 \>2 



or 



A 2 =p 2 -4p 3 > 4|ps|, (22) 

Since ^2 > is supposed, the condition lt%2l reads P3 < or p 2 > 8p 3 > 0. 
Thus, we obtain the following conditions for our solutions to be stable. 



P3 + PiPa = 

p 3 < or p 2 > 8p 3 > 0. 



(23) 



5 Numerical analysis of the system 

First, we we conduct numerical analysis of the system with random initial per- 
turbations in order to confirm validity of the stability condition (|23|l . We set 
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Figure 5: vs. if for stable Figure 6: a+(K) vs. if for unsta- 

parameters (pi = —0.5, P2 = 4.0, ble parameters (pi — 1.3, P2 = 3.0, 
p 3 = 0.5, ;p 4 = 1.0). p 3 = 0.25, P4 = 0.8). 

a boundary condition which enables our exact solutions to persist. It is sup- 
posed that, outside a boundary, there are oscillators whose phases are different 
from those just inside the boundary by w. This boundary condition is called a 
vanishing boundary condition in this paper. For values of parameters satisfying 
the stability condition 123(1 . both target patterns and spiral patterns are found 
to be stable even for random perturbations. For example, in a stable case that 
Pl = —0.5, p 2 = 4.0, p 3 — 0.5 and p 4 = 1.0, a single target pattern and the 
pattern with two spirals are stable even for max 20% random perturbations as 
shown in Fig[7| 
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Figure 7: Initial 20% random perturbations are added to exact solutions in a 
stable case that pi = — 0.5,^2 = 4.0, p$ = 0.5, p^ = 1.0. Length of an arrow 
indicates \W n \ and its direction means argW n . 



Next, we show that a single target pattern can be reached from some initial 
conditions far from the exact solution. For example, setting initial values as 
W n = 1/2, <j)(n) = 7r/4 for all n in the case pi = 0.5, P2 = 3.0, pz = 0.25 
and p4 = —0.5, we eventually have a target pattern, where \W n \ — 2, shown 
in FigEII We also find other initial conditions far from the exact solution that 
lead to the single target pattern but those do not give enough informations for 
a basin structure of the target patterns. 
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Figure 8: A target pattern of the numerical solution at time t = 50 for the 
initial condition given in Fig0and \W n \ 2 vs. t for all n, which shows stability 
of the target pattern. 




Figure 9: A combined spiral pattern of the numerical solution at time t = 50 
for the initial condition given in FigQand |Wn| 2 vs. t for all n, which shows 
stability of the combined spiral pattern. 




Figure 10: An initial condition \W\ = 0.5 and phase= ir/4 for all oscillators and 
a target pattern at time t — 150 developed from the initial condition. 
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6 Van der Pol's Oscillators 



Similar exact target and/or spiral solutions can be also constructed for other 
diffusively and dispersively coupled systems of oscillators which undergo the 
subcritical Hopf bifurcation. For example, we consider a system of van der 
Pol's oscillators which are coupled with each other diffusively and dispersively 
on a plane lattice. 

a 4 

Vn ~\~ ®X n ~Q*^ n ^ x ^ ; X n ) 
4 

-x n + D v y^(a?j - x n ), (24) 

where a is a positive constant and D x , D y are diffusive and dispersive constants 
respectively. Let us seek phase solutions of Ea. (|24f) similar to those discussed 
earlier. Assuming 

4 

X> = 0, (25) 
Eas. (|24l) are reduced to a single van der Pol's equation 

^X n + (Xl - l)j t X n +X n = 0, (26) 

where y/l + 4D y t -> t, x n = (1 - AD x /a) 1/2 X n and 

(l+4D y )>0, (l-4Z> x /a)>0 (27) 

are assumed. If the van der Pol's attractor of Ea. 1)26(1 is written as x n (t) — A(t), 
its opposite-phase function x n (t) = —A(t) is also the attractor of Ea. 1(26(1 . By 
choosing x n (t) as A(t) or —A(t) for all n so that x n (t) are satisfy the relation 
l(25|) , we can construct various exact solutions whose phase patterns are exactly 
the same as those in the discrete CGL equation. However, the stability analysis 
will be much more difficult than the previous case. 



dt 

d_ 

dt 



Vn 



7 Conclusions 

We present various exact solutions of a discrete CGL equation on a plane lattice, 
which describe target patterns and spiral patterns and derive their stability 
criteria. Exact solutions obtained here represent not only a single target pattern 
and a single spiral pattern but also patterns consisting of targets and spirals. In 
the center of a spiral pattern, at least one defect is necessary and a number of 
spiral's arms depends on a number of defects. For example, a defect produces 
a two-arm spiral and four defects in the center generate a four-arm spiral and 
so on. For a system of the van der Pol's oscillators, we can also construct 
various exact solutions whose phase patterns are exactly the same as those in 
the discrete CGL equation, pp 



10 



Acknowledgement 



One of the authors (K.N.) has been, in part, supported by a Grant-in-Aid for 
Scientific Research (C) 13640402 from the Japan Society for the Promotion of 
Science. 



References 

[1] K.Nozaki and N.Bekki, J. Phys. Soc. Jpn. 53 53 (1984), 1581. 

[2] I.S. Aranson, L. Aranson, L. Kramer, and A. Weber, Phys. Rev. A 46, 
R2992 (1992) 

[3] Igor S. Aranson and Lorenz Kramer, Rev. Mod. Phys. 74, 99 (2002) 

[4] J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems 
and Bifurcations of Vector Fields, Springer- Verlag (1983). 

[5] S. Mori: M. Thesis, Faculty of Science, Ochanomizu University, Otsuka, 
Bunkyo-ku, Tokyo, 1996. 

[6] H. Yamada and T. Nakagaki, talk at Traffic and Granular Flow '01, October 
15-71, 2001, Nagoya University, Japan. 

[7] S. H. Strogatz, Nature 410 (2001), 268. 



11 



